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ABSTRACT 



The formation rate of luminous galaxies seems to be roughly constant from 
z ~ 2 to ~ 4 from the recent observations of Lyman break galaxies (LBGs) 
(Steidel et al. 1999). The abundance of luminous quasars, on the other hand, 
appears to drop off by a factor of more than twenty from z ~ 2 to z ~ 5 (Warren, 
Hewett, & Osmer 1994; Schmidt, Schneider, & Gunn 1995). The difference 
in evolution between these two classes of objects in the overlapping, observed 
redshift range, z = 2 — 4, can be explained naturally, if we assume that quasar 
activity is triggered by mergers of luminous LBGs and one quasar lifetime is 
~ 10 7 ~ 8 yrs. If this merger scenario holds at higher redshift, for the evolutions 
of these two classes of objects to be consistent at z > 4, the formation rate of 
luminous LBGs is expected to drop off at least as rapidly as exp (—(-2 — 4) 6 / 5 ) 
at z > 4. 



Subject headings: Cosmology: large-scale structure of Universe - cosmology: 
theory - quasars 
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1. Introduction 



Observations of galaxies in the rest frame UV band ( [Lilly et al. 1996; Madau e~ 



|al~|1996; |Connolly et~al~| 1997; [Sawicki, Lin, fc Yee| 1997; [lreyer et aL~| 1998; 



Pascarelle, Lanzetta, fc Fernandez-Soto| 1998) indicate that the galaxy formation 



rate rises steeply from z — to z ~ 1, with a nearly constant rate thereafter up to z ~ 4 
( ISteidel et aL~| l999). While at low redshift (z < 2) the evolution of luminous quasar 
abundance resembles that of luminous galaxies (e.g., Sanders & Mirabel 1996; Boyle & 
Terlevich 1998), at high redshift [z > 2) the two classes of objects do not seem to parallel 
one another, with the luminous quasar formation rate (e.g., Warren et al. 1994; Schmidt et 
al. 1995) dropping off more steeply than that of luminous galaxies. 

In this Letter a phenomenological approach is taken to relate the observed formation 
rate of luminous LBGs to the observed abundance evolution of luminous quasars at 
z > 2. It is shown that, if 1) quasar activity is triggered by LBG mergers and 2) each 
quasar period lasts ~ 10 7 ~ 8 yrs, then the apparent difference, in both shape and amplitude, 
between the evolutions of bright LBGs and bright quasars from z = 2 to z = 4 can be 
explained quantitatively. The first assumption finds its support from both the observational 
evidence that a significant fraction of quasar hosts have disturbed morphologies or ongoing 
galaxy-galaxy interactions (e.g., Boyce et al. 1996; Bahcall et al. 1997; Boyce, Disney, & 
Bleaken 1999) and the theoretical consideration that merger of two (spiral) galaxies seems 
to provide a natural mechanism to fuel the central black hole (e.g., Barnes & Hernquist 
1991). The second assumption is also theoretically well motivated (Rees 1984, 1990) and 
now strongly implied or required by the mounting observational evidence that most nearby 
massive galaxies seem to harbor inactive black holes at their centers (e.g., Richstone et al. 
1998). 

The primary purpose of this work is to use this merger model to infer the LBG 
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formation rate at higher redshift z > 4. Given the precipitous drop-off of luminous quasar 
(M B < —26.0) abundance from z = 2 to z = 5 the formation rate of luminous [Mab > —23 
to —22) LBGs at higher redshift (z > 4) is predicted to drop off as least as fast as 
oc exp 

4 ) 6 / 5 ), if the merger scenario holds. A cosmological model with go — 0.5 and 
Hubble constant H = 50km/sec/Mpc is assumed for the analysis presented here. 

It is noted that this simple merger model would probably fail at z < 2 without having 
taken into account the evolution of gaseous fuel supply to the central black holes in galaxies 
(Kauffmann & Haehnelt 1999; Haiman & Menou 1999). 



2. Galaxy Merger Rate and Quasar Abundance 

Denoting f(z) as the galaxy formation rate (galaxy formation per unit time per unit 
co moving volume) as a function of redshift, then the (cumulative) number density of formed 
galaxies (number of galaxies per unit comoving volume) is (ignoring the small fraction of 
galaxies that merge) 

9w=rjw^. (i) 

For simplicity Qo = 1 will be assumed, which should be a good approximation at high 
redshift (z > 2) for the range of cosmological models of current interest (fl > 0.2). The 
merger rate for a galaxy with an internal one-dimensional velocity dispersion (Ji{z) in a 
cluster/group with galaxy number density d(z) (assuming all galaxies under consideration 
are identical) and one-dimensional velocity dispersion a e (z) has been computed by Makino 
& Hut (1997, equation 33) to be: 

18 1 

P(z) = ^—d(z)r v (zfa t (z)R(x), (2) 

where R(x) is a dimensionless function of x(z) = o e jui which depends on the galaxy model 
and r v (z) is the virial radius of a galaxy. Makino & Hut (1997) demonstrate that R(x) is a 
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constant (~ 11 — 14) to good accuracy for x > 2 for several different galaxy models. 

Clearly, not all galaxies participate in merging at any given time; most galaxies have 
merger time scales much longer than the Hubble time. Rather, only galaxies in dense 
environments such as groups or clusters of galaxies have significant probability to merge 
with others. To make the problem more tractable it is assumed that a fraction, fl(z), of all 
galaxies [g(z)] under consideration at any given time is in dense environments (i.e., typical 
groups/clusters at z) where most mergers occur, and the remainder of galaxies (i.e., field 
galaxies) have zero probability of merger. Then, the total merger rate is 



where £q(z) is the assumed quasar lifetime (assuming that tq is much less than the Hubble 
time, which turns out to be necessary for the model to be viable). 

There are two significantly uncertain remaining parameters, d(z) and a e (z), which 
need to be specified. It is noted that quasar activities at high redshift seem to occur mostly 
in regions with galaxy number density typical of present-day clusters/groups of galaxies. 
This information is provided by observations of quasar companions which have a typical 
separation from a quasar of a few hundred comoving kiloparsecs (e.g., Djorgovski 1999). At 
redshift z = 1 — 2 there is evidence from larger observational data sets that quasars reside 
in cluster-like environment (Hall & Green 1998). It thus appears that d(z) may be a weak 
function of redshift and is assumed to be constant here (more discussion on this later). The 
velocity dispersion of characteristic systems (groups/clusters in this case), u e (z), should 
be a decreasing function of redshift in any hierarchical cosmological model. Here we take 




(3) 



and the quasar abundance at any given redshift z is 




(4) 
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advantage of the insight of Kaiser (1986) and use the solution for simple power-law models: 

1 n— 1 

a e (z)=a e (0)(l + z)^, (5) 

where n is the power index of the primordial density fluctuation spectrum at the relevant 
scales for clusters/groups. For cold dark matter like models or from observations of local 
large scale structure n is expected to be ~ — 1. 

The purpose of estimating LBG formation rate at z > 4 is met by finding f(z) at z > 4 
that matches the observed quasar abundance in the range z > 2. For the present analysis a 
simple functional form of LBG formation rate is adopted: 

f(z) = A for 2 < z < 4 

f(z) = Aexp(-(^-4) 6/5 ) for z > 4, (6) 

consistent with the latest LBG observations at high redshift up to z = 4 (Steidel et al. 1999), 
where A is a normalization constant. At z > 4, where observations are unavailable, a simple 
form is proposed so as to provide an adequate fit to the observed quasar abundance at z > 4 
(see Figure 1 below). Using equations (1-3,5-6), we find Q(z) (equation 4), shown as the 
heavy solid curve in Figure 1. Here, for the shown Q(z) we use n = —1.0, cr e (0) = 10 3 km/s, 
<7j = 100km/s, P = 0.025 (being constant which is consistent with the adoption of n — —1 
powerlaw model), d = 40.0 /i 3 Mpc -3 , r v = 200/?," 1 kpc, R(x) = 12 and t Q = 3 x 10 7 yrs- A 
cosmological model with q = 0.5 and Hubble constant H = 50km/sec/Mpc is assumed. 
Also shown as symbols are observational data of bright quasars (M B < —26.0): open 
circles are from Warren et al. (1994) and solid dots are from Schmidt et al. (1995). The 
open square from Kennefick, Djordovski, & de Carvalho (1995) for M B < —26.7 quasars is 
shown to indicate the steepness of quasar luminosity function near the absolute magnitude 
M B ~ -26.0. 

It is seen that the merger model provides an adequate fit to the observed luminous 
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quasar abundance in the entire redshift range considered (z > 2). The dashed curve in 
Figure 1 shows f(z) with arbitrary vertical units. The dotted curve in Figure 1 shows g(z), 
normalized to be 1.0 x 10- 4 /i 3 Mpc" 3 at z ~ 3. Note that Figure 5 of Steidel et al. (1999) 
shows the differential luminosity function of UV bright LBG galaxies (i.e., star- forming 
galaxies), calling it gd(z), while here, g(z) is the cumulative density of formed galaxies. 
Roughly speaking, if f(z) is constant, then g(z)/gd(z) = tfj(z) /tsF, where tu{z) is the 
Hubble time at redshift z and tsF is the star (burst) formation duration (i.e., LBG phase) 
of a galaxy. Since t H (z)/t SF 10 s > yrs/10 8 yrs ss 10, the above normalization roughly 
corresponds to LBGs with gd(z) ~ 10~ 5 /i 3 Mpc~ 3 , which in turn corresponds to LBGs with 
M AB = -23 to -22 (Figure 5 of Steidel et al. 1999). 

3. Discussion 

On one hand, as Q(z) at z < 4 does not depend sensitively on the form of f(z) at z > 4, 
the good agreement between Q(z) and the observed quasar abundance in the redshift range 
z = 2 — 4 (where both types of objects are observed) suggests that the merger scenario of 
luminous LBGs provides a quantitatively viable model for bright quasar formation. On the 
other hand, Q(z) at z > 4 does depend sensitively on the adopted form of f(z) at z > 4. 
The fact that the proposed model yields an overall shape at z = 2 — 5 that fits observations 
implies that the luminous LBG formation rate should drop off at z > 4 as indicated by 
f(z) in eq. 1, if merger scenario holds at z > 4. But to have a secure estimate of f{z) 
at z > 4, it is vital to understand the dependences of Q(z) on various other parameters, 

3 n — 1 

namely, Q(z) <x (3(z)d(z)af(z)r%(z)tQ(z)(l + z)~z~^+z. We have set each of the parameters 
constant (independent of redshift), which is considered to conservative in the following 
discussions if a more likely redshift dependence of the quoted parameter (holding all other 
parameters constant) would require an even steeper decreasing function for f{z) at z > 4 
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than indicated by equation (6). Let us now examine each parameter to assess how each 
parameter may vary with redshift. 

First, it seems that cn(z), r v (z) and tQ are likely to decrease with redshift, making the 
assumption of their being constant conservative. 

Second, (3 = 0.025 is equivalent to the assumption of mergers taking place in galaxy 
systems corresponding roughly to 2a peaks and has implications for the correlation function 
of quasars. The bias factor of halos over mass is b = 1 + (y 2 — 1)/S C (Mo & White 1996), 
equal to 2.91 for v = 2 and 5 C = 1.57. If the cluster-cluster correlation function has a 
shape oc r~ 2 (close to the usual slope of —1.8), then the correlation length of clusters is 
br m , where r m is the correlation length of the underlying mass and evolves as oc (1 + z)^ 1 
(Kaiser 1986) for n = —1 and Q = 1- Our choice of f3 — 0.025 consequently implies a 
correlation length for quasars of approximately 2.91r m (0)/(l + z), which is equal to ~ 5/i _1 
comoving Mpc at z ~ 2 (using r m (0) ~ 5.0/i _1 Mpc), in agreement with what is observed for 
quasars (e.g., Kundic 1997; Boyle et al. 1998). In any case, it is unlikely that (3 decreases 
with redshift. Therefore, setting (3(z) constant is conservative. An important implication 
of this model is that the comoving correlation length of luminous quasars should decrease 
with redshift no faster than (1 + z)^ 1 at z > 2, a potentially testable prediction. Stephens 
et al. (1997) give a correlation length of z > 2.7 quasars of 17.5 ± 7.5/i~ 1 Mpc. It will be 
very valuable to determine the correlation length of high redshift quasars with significantly 
smaller errorbars. 

Third, observations may have indicated that d(z) may be an increasing function of 
redshift at z > 4 (Djorgovski et al. 1997; Djorgovski 1999). Therefore, assuming d(z) to be 
constant is conservative. 

Finally, for a plausible power spectrum (such as CDM like) n is likely to be smaller 
at smaller scales thus smaller at higher redshift. Thus, assuming n to be a constant 
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is conservative. Overall, our assumption of constancy for various parameters seems 
conservative; i.e., f(z) should decrease at least as rapidly as indicated by equation (6) at 
z > 4. 

All the analyses so far have been based on the available (optical) observations of 
quasars, which appears to indicate a sharp drop-off of quasar abundance at z > 4. Dust 
obscuration effects are often invoked to explain the apparent drop-off of quasar abundances 
at high redshift (e.g., Ostriker & Heiler 1984; Pei 1995). However, recent radio surveys of 
high redshift quasars seem to indicate that the drop-off of the number density of bright 
radio quasars is very similar to that from optical surveys (e.g., Hook, Shaver, & McMahon 
1998) with the implication that the effect of dust on the observed drop-off of bright quasars 
at z > 2 may be small. 

One potential problem with the merger model is that observations show that a large 
fraction of quasar hosts at low redshift (z < 0.5) appear to be quite normal looking, 
i.e., without disturbed appearances. But one would expect that, if galaxy-galaxy merger 
time scale is longer than the proposed quasar lifetime, all quasar hosts should display 
appearances of some interaction. One possible solution to this problem is that quasar 
formation is delayed, i.e., a quasar does not start to shine until the galaxy-galaxy merger 
is nearly complete. In other words, the time it takes to set up the central (BH) region for 
quasar activity during galaxy merger may be comparable to the time that it takes for the 
two galaxies to merger. 

4. Conclusions 

In an early classic paper Efstathiou & Rees (1988) show that quasar abundance at high 
redshift can be accounted for in the standard cold dark matter model if massive halos are 
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related to the formation of black holes, with an intriguing prediction that the abundance 
of luminous quasars should decrease rapidly beyond z = 5 (for a more recent treatment 
see Haehnelt & Rees 1993). [The evolution of low-luminosity quasars, of course, does not 
necessarily have to follow that of their luminous counterparts (e.g., Haiman & Loeb 1998)]. 

In this Letter a different approach is taken by directly relating the observed evolution 
of luminous LBGs to the observed evolution of luminous quasars at high redshift (z > 2). 
With a set of seemingly reasonable parameter values, it is shown that consistency between 
the two classes of objects at z = 2 — 4, where both classes are observed, can be achieved, if 
one assumes that 1) Lyman break galaxies merger to trigger quasar activity and 2) quasar 
lifetime is ~ lO 7 " ~ 8 yrs. At z > 4, consistency can be achieved, only if additionally the 
formation rate of luminous LBGs drops off as exp(— [z — 4) 6 / 5 ) or faster, a prediction that 
may be tested by future observations. One implication from this model is that LBGs with 
Mab > —23 to —22 merge to form quasars with M B < —26.0 at z > 2. Correlation analysis 
of relevant LBGs and quasars should shed light on this. 

At lower redshift additional, more model dependent assumptions regarding the 
supply of available gas to fuel black holes would be required to make qualitatively viable 
predictions. Kauffmann & Haehnelt (1999; see also Haiman & Menou 1999) have presented 
a detailed model, based also on merger scenario, to unify the evolution of galaxies and 
quasars in the cold dark matter model under several plausible assumptions concerning the 
evolution of fuel gas to the central black holes. The success of the model of Kauffmann 
& Haehnelt (1999) at low redshift [z < 2) and the model presented here at high redshift 
(z > 2), both based on galaxy merger scenario, suggests that galaxy merger may play an 
indispensable role in quasar formation. 
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Fig. I. — The heavy solid curve shows the evolution of the abundance of bright quasars 
computed using the merger model of this paper. Also shown as symbols are observational 
data of bright quasars {Mb < —26.0): open circles are from Warren et al. (1994) and 
solid dots are from Schmidt et al. (1995). The open square from Kennefick, Djordovski, & 
de Carvalho (1995) for M B < —26.7 quasars is shown to indicate the steepness of quasar 
luminosity function near the absolute magnitude Mb ~ —26.0. The dotted and dashed 
curves show g(z) (eq. 1) and f(z) (eq. 6), respectively. 
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